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In this paper we study the problem of the existence of a least-action principle for 
invertible, second-order dynamical systems, discrete in time and space. We show that, 
when the configuration space is finite, a least-action principle does not exist for such 
systems. We dichotomize discrete dynamical systems with infinite configuration spaces 
into those of finite type for which this theorem continues to hold, and those not of finite 
type for which it is possible to construct a least-action principle. We also show how 
to recover an action by restriction of the phase space of certain second-order discrete 
I dynamical systems. We provide numerous examples to illustrate each of these results. 
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1 Introduction 



^ . Hamilton's Principle of Least Action provides a unifying and elegant reformulation of classical 
mechanics. If we suppose that the instantaneous state of a mechanical system is described by 
generalized coordinates r taking values in a configuration space C, then Hamilton's Principle 
^ ' states that the trajectory of the system between times and T will be such as to extremize 
<^ • the action functional 

X: S[r]= I dtC{r,r), (1) 







where the Lagrangian C : TC — > M is a real-valued function on the tangent bundle of C, and 
where the dot denotes differentiation with respect to time t. Setting the Frechet derivative 
of S[r] to zero, 

o=^^=^-if^y (2) 



5r{t) dr dt \dr J 
yields Lagrange's equations of motion for the mechanical system. 

Example 1 A mechanical system with mass matrix M and potential energy function V : 
C — s> M has the Lagrangian function 

C{r,r) = ^r^Mr -V{r), (3) 
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where the T superscript denotes transpose. Lagrange's equations of motion for this system 
are easily seen to be equivalent to Newton's equations 

dV 

Mr = - — . 
or 

Newton's equations are second-order autonomous ordinary differential equations that are 
well known to be reversible; that is, they are invariant under time reversal, t — > —t. They 
can be reformulated as a system of first-order differential equations known as Hamilton's 
canonical equations. 



dH 
dp 

dH dV 



M-'p (4) 



^ dr dr' ^ ^ 



where we have introduced the Hamiltonian H = ^p'^M^^p + V{r) and the momentum p. 
Under time reversal the momentum changes sign, p —p, and the time derivatives change 
sign, d/dt — > —d/dt, so Hamilton's canonical equations are invariant. 

All of the above is well known for mechanical systems with continuous configuration 
spaces C, and with a continuous time coordinate t. If time is discrete, as it is in a finite- 
difference algorithm, the configuration of the system at time t G Z may be denoted by r{t), 
and the trajectory of the system from time to time T G Z+ is given by the sequence 
r = {r(t) I < t < T}. In this case, we might suppose that the action functional in Eq. ([T]) 
could be recast as a function of r, 

T-l 

S{r)=J2C{r{t),r{t + l)), (6) 
t=o 

where £ : — > M is the discrete-time analog of the Lagrangian function. The equations of 
motion may be obtained by setting to zero the gradient of the above, 

0=1^= A(r(t-l),r(t)) + A(r(t),r(t + l)), (7) 

where £i and £2 denote the gradient of C with respect to its first and second arguments, 
respectively. If £ is a symmetric function of its two arguments, then \/a,b£C '■ 'C.i{o,,b) = 
C2{b, a), from which it quickly follows that Eq. (JTj) is invariant if t — 1 and t + 1 are inter- 
changed, demonstrating reversibility of the discrete-time mechanical system. This formula- 
tion of discrete-time mechanics seems to have been first posited by Baez and Gilliam [1]. 
More recently, it has recently been demonstrated [2] that this approach may be used to 
construct symplectic finite-difference algorithms for the numerical simulation of mechanical 
systems. 

The utility of action principles for mechanical systems in discrete time with a finite 
discrete set of states, such as cellular automata, has received much less attention. We might 
still suppose that Eq. ([6]) is a reasonable definition for the action, but it must be understood 
that r(t) takes its values in a discrete configuration space C of finite cardinality, such as a 



2 



finite set of integers. One immediately runs into tlie difficulty that it is impossible to take 
derivatives of £ with respect to its two arguments, so £i and £2 are not well defined. Baez 
and Gilliam addressed this problem using methods of algebraic geometry [1], but here we 
take a different approach. Rather than attempt to salvage Eq. ([71), we instead give primacy 
to Eq. ([6]). For a discrete configuration space C, we pose the question of whether or not it 
is possible to find a real- valued function C : C"^ —>■ such that the action S{r) given by 
Eq. ([6]) is smallest for trajectories governed by a reversible evolution rule for the finite-state 
system- one giving r{t + 1) in terms of r{t) and r{t — 1) that is invariant under interchange 
of if: -|- 1 and t — 1. In this paper, we answer this question in the negative. Indeed, we 
disprove the existence of an action principle in a finite configuration space not only for the 
case of reversible dynamical systems, but also for the much more general case of invertible 
dynamical systems. The proof is general enough to apply to any invertible dynamical system 
with a finite numbers of states, including all invertible cellular automata. 

Our approach is based on an observation of Feynman [3] who noted that, if an action is 
minimized on a physical path, then it must be minimized on every "infinitesimal subsection" 
of that path. Because our system evolves in discrete time steps, we replace the notion of 
"infinitesimal subsection" with its discrete counterpart, examining the existence of action 
principles for a sequence of three time steps. This approach leads to a set of inequalities 
that must be satisfied by the Lagrangian function, and we show that these inequalities are 
necessarily inconsistent for any invertible dynamical system with a finite number of states. 

Because this paper is meant to be self-contained, we begin by reviewing the notions of 
reversibility and invertibility for discrete dynamical systems. We then describe a represen- 
tation of such systems via elements of the symmetric group; this representation will play an 
important role in our subsequent proof. Finally, we make precise the concept of action for a 
finite-state dynamical system and we present a proof of the "No-Go Theorem" for an action 
principle. 

We conclude the paper by observing that, if we relax the hypothesis of finiteness of the 
configuration space or if we work on a suitable restriction of the phase space, an action 
principle can be recovered. Two examples describing how an action can be recovered in 
these cases are discussed in the conclusions section of our work. 

2 Discrete Dynamical Systems 

We consider dynamical systems that are discrete in time and discrete in state. Discreteness in 
time means that the evolution of the system takes place in successive time steps. Discreteness 
in state means that at each time step the system may be in only one of r G Z"*" states; we 
refer to the set of all possible instantaneous states of the system as the configuration space, 
C = {0, 1, . . . , r — 1}. (Later in the paper, we consider the limit of infinite configuration 
space, r — > 00.) The state of the system at time t will be denoted by r(t). 

Example 2 The instantaneous state of many discrete dynamical systems, such as cellular 
automata, is most conveniently described by n dependent variables, each of which takes its 
values in a set Q, with \ Q\ = q. This is equivalent to describing the system by a single state 
variable which takes values in a configuration space C = Q" of size r = q^ . In this case, we 



3 



denote the j 'th dependent variable at time step t by rj(t) G C, for < j < n and t > 0, so 
that r{t) may be thought of as an n-vector. 

We refer to a dynamical system as p 'th order if the state at time t depends on those at 
times t — l,t — 2, . . . ,t — p. In this paper, we shall focus on the case p = 2. A first-order 
{p = 1) djTiamical system has an evolution rule of the form 

r(t+l) = /(r(t)), (8) 

where we have introduced the function f : C ^ C. A second-order [p = 2) dynamical system 
has an evolution rule of the form 

r{t + l)=g{r{t),r{t-l)). (9) 

where we have introduced the function g : —>■ C. More generally, a p'th-order dynamical 
system requires the introduction of a map — > C A p'th-order dynamical system has a 
well-posed initial-value problem if r is specified at p successive initial times, say from t = 
to t = p — 1. In what follows, we denote by F the dynamical system induced by g. 

2.1 Configuration Space, Phase Space, and Periodic Orbits 

For deterministic first-order dynamics, p = 1, it is evident that the system will return to 
an initial configuration in at most |C| = r time steps, after which the motion repeats. For 
p > 1, however, the evolution in configuration space is not manifestly Markovian, since the 
configuration of the system at time t + 1 depends on that at p earlier time steps. 

To recast the dynamics as manifestly Markovian, it is necessary to consider the set of all 
p-tuples of configurations. Each such p-tuple shall be called a phase point, and the set of all 
of them will be called the phase space, V = C^. The advantage of this approach is that the 
p-tuple at time t + 1 depends only on that at time t; that is, the phase-space dynamics are 
first-order. 

The system will return to an initial phase point in at most \V\ = \C\^ = time steps, 
after which point the motion repeats. The number of time steps required to return to an 
initial phase point is called the recurrence time and denoted by r. Generally, the recurrence 
time will depend on both the dynamical rule and the initial condition 

2.2 Invert ibility and Reversibility 

A p'th-order dynamical system is said to be invertible if it is possible to uniquely determine 
r{t — p) given r{t — p + l),...,r{t). For a first-order dynamical system, invertibility means 
that it must be possible to solve Eq. ([H]) for r(t) or, in other words, that / must be an 
invertible function. For a second-order dynamical system, this means that it must be possible 
to solve Eq. ([9]) for r{t — 1). Thus, an invertible dynamical system admits a rule for evolving 
backwards in time, as well as forwards in time. Alternatively stated, it has a well posed 
final- value problem, if r is specified at p successive final times, say from t = T — p+1 to 
t = T. 
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An invertible dynamical system is said to be reversible if the rule for evolving backwards 
in time is identical in form to that for evolving forwards in time [5j. This is substantially 
more restrictive than invertibility. For the first-order [p = 1) dynamical system in Eq. ([8]), 
it means that 

r(t) = /(r(t + l)). (10) 
For the second-order {p = 2) dynamical system in Eq. ([9]), it means that 

r{t-l) = g{r{t),r{t + l)). (11) 

These requirements place severe restrictions on the form of the functions / and g. For 
example, in the first-order {p = 1) case, we can combine Eqs. (IHl) and (ITOl) to find 

r{t) = f{f{rm. (12) 

Now any possible r(t) G C could be specified, perhaps as an initial condition, and the above 
must hold for any of them, so we require 

\/xec ■■ X = f if (x)) . (13) 

Thus, any reversible, first-order, discrete-time dynamical system must have a recurrence time 
of at most two. This makes them substantially less interesting than reversible dynamical 
systems of higher order. 

In the second-order (p = 2) case, we can combine Eqs. ([9]) and (ITT]) to find 

r{t-l)=g{r{t),g{r{t),r{t-l))). (14) 

Again, any possible 2-tuple {r(t — 1), r{t)} E V = could be specified, perhaps as an initial 
condition, and the above must hold for any of them, so we require 

yx,y&c ■■ X = g{y,g{y,x)) . (15) 

The requirements expressed by Eqs. (1131) and ( |T5|) are not as different as they might first 
appear, since y may be regarded as a parameter in the second equation; that is, we may 
write g{y,x) = gy{x). So in either case, we seek functions - either f or gy mapping 
C — s> C whose self-composition is the identity. If we find more than one such function for a 
given configuration space, we may let the first argument of g, namely y, determine which 
one we use for the second-order case. 

2.3 Permutation Representation 

Let F be a second-order discrete dynamical system. From the discussion above we have that 
F is invertible if the associated function g{y, x) = gy{x) is a bijection gy : C ^ C; that is, F 
is invertible if gy acts as a permutation on C for all y EC. We also have that F is reversible 
if gy is a bijection and self-inverse; that is F is reversible if g\^ is the identity, or gy acts as 
a permutation of order at most two on C for sX\ y E C so that it is a product of disjoint 
transpositions. The most obvious way to describe g is by explicitly listing the functions gy 
for all y E C. 
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Example 3 A dynamical system of order p = 2 with configuration space C = {0,1,2} of 
cardinality r = 3 may have dynamics described by the function 
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y 




9y 
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U 
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2 
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2 
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1 


e 


2 


2 


2 





T/ie last column of this table lists the elements of the symmetric group S3 used for each value 
of y; in this case, cxi = (021) is used for y = 0, and the identity permutation e is used for 
y = I and y = 2. Therefore the function g is uniquely specified by the sequence {ai,e,e}. 
Because ai has order three, this system is invertible but not reversible. 

Example 4 We consider the dynamical system of Example but this time with dynamics 
described by the function 
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y 


9vix) 


9y 
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To 


2 
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2 


1 


2 
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e 


2 


2 


2 





This example differs from the previous one in that tq = (21) is used for y = 0. Therefore 
this function g is specified by the sequence {ro,e, e}. Because tq has order two, this system 
is reversible. 

From the above examples, we see that we can always represent g uniquely as a sequence 
of r elements of the symmetric group Sr, each of them corresponding to a permutation gr^, 

9 = {go, 91, ■ ■ ■,gr-i} , 

where gj is now used to denote general elements of the permutation group Sr- If all elements 
of this sequence are identical, so that gy does not depend on y, we say that the associated 
dynamics F are trivial; otherwise, they are nontrivial. 
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3 Action Principle 

Definition 1 A path of length T connecting states and Vf is a sequence of states in 
configuration space, r = {r(t) G C | < t < T}, beginning with state r{0) = Vi and ending 
with state r{T) = Tf. 

The following definitions refer to second-order dynamics F described by a function g: 

Definition 2 We say that a path connecting r(0) = and r(T) = rj is a physical path 
of length T if either T G {1, 2}, or T > 2 and r{t) = g{r{t - 1), r{t - 2)) for 2 < t < T . 
Otherwise, r is said to he a non-physical path. 

That is, a path is physical if every sequence of three consecutive points on it, {r{t — 2), r{t — 
l),r(t)} obeys the evolution equation, Eq. ([9]). 

Given any pair of states Vi^Vf G C, let us further denote the set of all physical paths 
between and Vf by f/(rj; rj), the set of all the non-physical paths between and rj by 
V{ri] Vf), and their union by W{ri; Vf) = U{ri] rj) U V{ri; rj). 

Definition 3 Given a Lagrangian function C : CxC — > M, the action of a path r is defined 
by 

T-l 

S{r) = Y,^r{t),r{t+l)). (16) 
t=o 

Definition 4 We say that T admits an action principle if, for any ri,rf G C such that 
U{ri, Vf) 7^ 0, and for any r G U{ri; Vf) and r' G V{ri; Vf), it is true that 

S{r) < S{r') (17) 

for some Lagrangian function C : C x C — > M. 

That is, the action on any physical path is required to be less than the action on any 
nonphysical path. We do not require any relative ordering of the action on two physical 
paths. 

The following Lemma will be central to our analysis: 

Lemma 1 Let r G (/{rf, r f) he a physical path of length T . IfV admits an action principle, 
then 

C{r{t), r{t + 1)) + C{r{t + 1), r{t + 2)) < C{r{t), r'{t + 1)) + C{r'{t + 1), r{t + 2)) (18) 
for all r' (t + 1) E C such that 

gir'{t + l),r{t))^gir{t + l),r{t)). (19) 

Proof. We first notice that 

r' = {r, = r(0), r(l), . . . , r(t), r'{t + 1), r(t + 2) . . . , r(T) = r^} 
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is a non-physical path of length T. Since, by hypothesis, we know that S{r) — S{r') < 0, we 
have that 

S{r) - S{r') = 

L{r{t), r{t + 1)) + C{r{t + 1), r{t + 2)) - C{r{t), r\t + 1)) - C{r\t + l),r{t + 2)) < 0. 



which proves the claim. 



□ 



Theorem 1 Invertible, second-order, nontrivial, discrete dynamical systems with a finite 
configuration space do not admit an action principle. 

Proof. Let F be a second-order discrete dynamical system with associated function g repre- 
sented by the sequence of permutations {fi'ojfi'i? • • • jfi'r-i}- Because the dynamical system is 
nontrivial, we may suppose that there exist i, j £ C such that Qi ^ Qj. If we set p = gjg^^ ^ e, 
we can write gj in terms of gi as 

9j = P9i- 

Using the tabular representation developed in the previous section, we can describe and 
gj as follows: 



X y 


9i{x) 


9y 


i 






1 i 






2 i 




9i 


r — 1 i 


9i{r - 1) 





X y 


9iix) 


9y 
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p{9^m 




1 J 


Pi9^{l)) 




2 j 


p{9m) 


9j 


r-1 j 


P{9iir-l)) 





In order to simplify the notation, we also define = 9i{k). 

To illustrate the main idea behind the proof of the theorem, we first consider the case 
when p is a transposition; that is, for some h,k & C, we have that p = (hk) (in cycle 
notation). Then the tabular representation of g^ and g^ can be written as follows: 
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9iix) 


9y 





i 
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i 
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a 


9i 
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i 






r — 1 
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9iix) 


9y 
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9j 
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a 




r — 1 


J 
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By analyzing the representation of g^, we notice that {h,i,^h} and {k,i,^k} are physical 
paths whereas {h,j,^h} and {k,j,^k} are not. Similarly, by looking at the representation of 
Qj, we notice that {h,j,^k} and {k,j,C,h} are physical paths whereas {h,i,^k} and {k,i,C^h} 
are not. Therefore, if F admitted an action principle, then by Lemma ([T]), C would satisfy 
the linear constraints 

C{h,z) + C{i,^f,) < C{h,j) + C{j,^f,) 
C{k,t) + C{t,^k) < C{k,j) + C{j,^k) 

c{h,j) + c{j,^k) < c{h,t) + c{t,^k) ■ 

C{k,j) + C{j,^h) < C{k,z) + C{t,^h) 

Adding the above inequalities, we notice that all terms cancel. Since the inequalities are 
strict, we arrive at a contradiction. 

Now let us generalize to the case where p E Sn- The tabular representation of Qj^ and Qj 
may be written as follows: 



(20) 
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9i{x) 


9y 


i 






1 i 






2 i 




9i 


r — 1 i 







X y 


9iix) 


9y 
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P(6) 




1 J 






2 J 


P(6) 




r -I j 


Pi^r-l) 





If we define Ip = {i E C \ p(^i) 7^ ^i}, then, for any k G /p, by Lemma [U we obtain the 
following contraints for the Lagrangian function C: 



C{k,i) + C{i,ik)<C{k,j) + C{],ik) 
C{k, j) + C{j, p(a)) < C{k, z) + C{t, p(a)) ■ 

If we add both sides of the fl2T|) . after cancellations we get: 



(21) 



f22) 



where k E Ip. 

After noticing that 



5^£(p(^,a)) = $^/:(^,efe) 



kelp 



kelp 



for any i E C, we can conclude that the sum of the left-hand sides of the inequalities 
coincides with the sum of their right-hand sides. Again, since the inequalities are strict, we 
arrive at a contradiction. 
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Example 5 Let us consider the discrete dynamical system T described by the following as- 
sociated function g : 
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// r admitted a non-trivial action principle, then, by looking at Qq and g^, we would infer 
that there exists a Lagrangian C verifying the following inequalities: 



£(0,0 
£{1,0 
£(2,0 
£(3,0 
£(0,1 
£(1,1 
£(2,1 
£(3,1 



+ £(0,1 
+ £(0,0 
+ £(0,3 
+ £(0,2 
+ £(1,2 
+ £(1,3 
+ £(1,0 
+ £(1,1 



< £(0,1 
<£(1,1 
<£(2,1 

< £(3,1 

< £(0,0 

< £(1,0 

< £(2,0 

< £(3,0 



+ £(Li) 
+ £(1,0) 
+ £(1,3) 
+ £(1,2) 
+ £(0,2) 
+ £(0,3) 
+ £(0,0) 
+ £(0,1) 



If we add the inequalities, after cancellations we arrive at the contradictory inequality < 0. 

4 Recovering an action 

We have proved that invertible second-order discrete dynamical systems with a finite con- 
figuration space C do not admit an action defined on the entire phase space — V. If we 
either allow C to be infinite, or if we restrict to a subset of V, the existence of an action 
principle is not ruled out by the above argument. In this section we present two examples 
of such instances where an action principle can be recovered. 

4.1 Infinite alphabet 

Using the notation previously introduced, let us consider a configuration space C of infinite 
size. Let then T be described by a function g that can be represented by a sequence of 
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automorphisms gi : Z — > Z 

9 = {gi}iez 

Definition 5 T is said to be of finite type if there exists i G Z+ such that gi contains some 
non trivial finite cycle. 

The following lemma is direct consequence of Theorem (1): 

Lemma 2 IfVis of finite type, then it does not admit an action principle. 

By the next example we would like to highlight the fact that, if F is not of finite type, 
an action principle may exist. 

Example 6 Consider the dynamics induced by g{y, x) = x — y, and let 

y) = ^ [{x + y)-{x + yf] + {x- yf 

with x,y G Z 

For any y' ^ y, we have that 

2 

C{x,y) + C{y,y -x) - [C{x,y') - C{y',y-x)] = --[-1 + 2{y - y')]{y - y'). 

We notice that the above expression is always negative since, if y — y' > then 

-l + 2{y-y')>0 
(since y — y' > 1). By the same token, if y — y < 0' , then 

-l + 2{y-y')<0. 

Hence we can conclude that C{x,y) is a Lagrangian for the dynamics defined by g{y,x) = 
x-y. 

The above action is, however, not unique. Indeed, let us assume a quadratic form for the 
Lagrangian C{x,y), that is 

C{x, y) = ax? + (5xy + 7?/^, 
with a, /9, 7 G M. For any y' 7^ y, we require that 

y)+C{y,y-x)- [C{x, y') - C{y\ y-x)] = {y- y') [y'{a + 7) + y(« + /3 + 7)] < 0. (23) 

It is easy to show that solutions to ff^gj) must obey following form: 

a + j3 = C a + /? + 7 = -C, 

where C (zM. is a constant. That is, C{x, y) = ax^ + (5xy + 7?/^ is an action for g for any 
a, /?, 7 G M such that 

2a + /? + 27 = 0. 
11 



4.2 Restricting the phase space 



Let r be a second-order discrete dynamical system, g its associated function, and let (x, y) G 
= "P be a point in its phase space. Then, as we have discussed in an earlier section, the 
dynamics in phase space V is represented by the map F: 

T -.V — >V 

{x,y) ^ {y,g{y,x)). 
We define the orbit originated at s G P as: 

a = {r\s) 1 1 G z+}, 

and decompose the global dynamics in phase space into a transient part 

W = {s eV \ y t>0 : T\s) ^ s} 

and a cyclic part 

Z = V\W 

Notice that, if F is invertible, then the transient part of the orbit through any point 
s G P is empty, and therefore W = ^. On the other hand, if F is non-invertible, then 
there are points in V whose orbits have a non-trivial transient part, and therefore W ^ ^. 
In this latter case, since the No-Go Theorem does not apply to non-invertible systems, an 
action may exist and we have explicitly computed a Lagrangian C in some of these cases 
(see Example [8] below) . 

If £ is a Lagrangian inducing an action for the non-invertible system F, then we can 
recover an 'invertible' action by restricting C to the truncated phase space Z = V\T. It is 
important to notice that this fact does not contradict our main result. In order to clarify 
that, given a finite set A, let us first define 

A* = {a:A^ — > A^ \ a{x,y) = {y,z), a G ^i^p} 

for any x,y & A and some z & A. Here denotes the symmetric group on elements. 

Example 7 Let A = {0,1}. Then = {00,01,10,11}, and it is not difficult to show 
that, out of the 4! possible permutations of elements of Al' , only four of them obey the 
above conditions, namely A* = {o"i}i=i,...,4, where ai{x,y) = {y,x), a2{x,y) = {y,l — x), 
(73 (x, y) = {y,x (By), 0-4(0;, y) = {y,l — x Q) y). Here x ®y denotes addition modulo 2. 

Now we can observe that, by the No-Go Theorem, F cannot be recast as a second-order 
dynamical system on a finite configuration space. That is, there exists no alphabet C such 
that F G C*. An interesting consequence of this observation is that, every time an action for 
some non-invertible system F is detected, we can recover some "negative information" on its 
limit cycle structure - namely, the fact that F cannot belong to C* . 

The following is an example of a non-invertible system which admits an action principle. 
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Example 8 Let \C\ — ?> and consider the second-order, non-reversible dynamics associated 
with the following function g: 



X 


y 


9v{x) 








2 


1 





1 


2 


C) 


(J 





1 


2 


1 


1 


1 


2 


1 








2 


1 


1 


2 





2 


2 






In phase space, the dynamics is given by a 7- cycle 

00 ^ 02 ^ 21 ^ 10 ^ 01 ^ 12 ^ 20, 

a 1-cycle 

11 ^ 11, 

and a single transient state 

22 ^ 20. 

We can compute the constraints for the Lagrangian with the help of the following di- 
agrams, which represent all possible paths (x(t),x(t+l),x(t+2)). More precisely, the left 
and the right corners of the diamonds represent states of the system at time t and at time 
t + 2., respectively. That is, they denote x{t) and x{t + 2), respectively. Then, for each pair 
{x{t) ., x{t + 2)) there are \C\ = 3 possible intermediate states x{t+ 1), which are represented 
by the column between x(t) and x(t + 2). Finally, we adopt the convention that plain arrows 
denote physical paths, while wiggly arrows denote unphysical paths. 

We begin with 








O-^l-^^O 2^^1^^2 2^^1^-1 i^^i^^2. 




2 2 2 2 



In the above diagrams, all paths are unphysical, hence these cases give rise to no constraints 
at all. 

Next, we consider 








2 
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In this case we notice that all paths are physical, therefore this case gives rise to no constraints 
as well. 

Next, we consider consider the cases where there is only one physical path connecting x{t) 
with x{t + 2); 








2 2 



The inequality constraints arising from the above cases are 

jC{0, 2) + jC{2, 1) < C{0, 1) + jC{1, 1) £(1, 2) + C{2, 0) < 1) + £(1, 0), 

£(0, 2) + £(2, 1) < £(0, 0) + £(0, 1), £(1, 2) + £(2, 0) < £(1, 0) + £(0, 0). 

Finally, we have to consider the cases when there are two physical paths connecting x{t) 
with x{t + 2): 






>1 ^2 1 *-l ^1 




2 2 



The inequality constraints arising from these cases are 

£(0, 0) + £(0, 2) < £(0, 2) + £(2, 2), £(1, 0) + £(0, 1) < £(1, 2) + £(2, 1), 

£(0, 1) + £(1, 2) < £(0, 2) + £(2, 2), £(1, 1) + £(1, 1) < £(1, 2) + £(2, 1). 

More generally, note that if we denote by Np the number of physical paths, and by Nn = 
\C\ — Np the number of unphysical paths, then the total number of inequalities is given by 
NpNn^Np{\C\-Np). 

As we have pointed out earlier, the No-Go Theorem does not preclude this non-reversible 
dynamics from having an action. In particular, if we remove the transient state 22 from 
the phase space, an action will be induced by a solution C{x,y) to the system of inequalities 
given by 

' £(0, 2) + £(2, 1) < £(0, 1) + £(1, 1) 

£(0, 2) + £(2, 1) < £(0, 0) + £(0, 1) 
£(1,2) + £(2,0)<£(1,1)+£(1,0) 
£(1, 2) + £(2, 0) < £(1, 0) + £(0, 0) 
£(1,0) + £(0,1) <£(1,2)+£(2,1) 
^ £(1,1) + £(1,1)<£(1,2)+£(2,1) 
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It can be seen that this system has an infinity of solutions ( therefore the action is not 
unique). A particular solution is given, for instance, by 

£(0, 0) = C{Q, 1) = £(0, 2) = C{1, 0) = £(1, 1) = 
£(1,2) =2 
£(2,0) = -3 
£(2,1) = -1. 

5 Conclusions 

In this work wc have discussed the question of the existence of action principles for invcrt- 
ible second-order discrete dynamical systems. We have proved that systems with a finite 
configuration space do not admit an action on the entire phase space. By weakening one of 
these two assumptions, that is assuming a configuration space of infinite size or working on a 
subset of the phase space, an action principle can be recovered and we have given examples 
of an action principle constructed in this way. 

Most of the motivations behind our work has come from the search for an action princi- 
ple for finite-state cellular automata. Given their discrete nature, we decided to avoid any 
reference to differential structures, and have provided a genuinely discrete definition of vari- 
ational principle by replacing the minimization of a continuous action with a system of linear 
inequalities obtained by minimizing a discrete action. In some sense, our result highlights 
both the similarities and the differences between continuous and discrete dynamical models. 
It seems to us that further investigations of these similarities and differences is warranted. 
After all, the elegance and power of continuous models in mechanics should not preclude the 
search for alternative, discrete models which could lead to useful applications, from both a 
theoretical and numerical-computational perspective. Future work might center on discrete 
analogs of other results in continuum mechanical systems, such as Noether's Theorem. 
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